Test of Jarzynski and Crooks fluctuation relations in an electronic system 



(N 

o 



(N 



c3 



i 

c 

o 
o 



(N 

>; 

OV 

O 
l> 

o 



X 



O.-P. Saira, 1 - 2 Y. Yoon, 1 T. Tanttu, 2 M. Mottonen, 2 ' 1 D. V. Averin, 3 and J. P. Pekola 1 

1 Low Temperature Laboratory, Aalto University, P.O. Box 15100, FI-00076 AALTO, Finland 
2 Department of Applied Physics/COMP, AALTO University, P.O. Box 14100, FI-00076 AALTO, Finland 
1 Department of Physics and Astronomy, Stony Brook University, SUNY, Stony Brook, NY 11794-3800, USA 

Recent progress on micro- and nanometer scale manipulation has opened the possibility to probe 
systems small enough that thermal fluctuations of energy and coordinate variables can be signif- 
icant compared with their mean behavior. We present an experimental study of nonequilibrium 
thermodynamics in a classical two-state system, namely a metallic single-electron box. We have 
measured with high statistical accuracy the distribution of dissipated energy as single electrons 
are transferred between the box electrodes. The obtained distributions obey Jarzynski and Crooks 
fluctuation relations within experimental error. A comprehensive microscopic theory exists for the 
system, enabling the experimental distributions to be reproduced without fitting parameters. 



Everyday concepts such as heat and mechanical work 
have their foundation in dynamics at the atomic and 
molecular scale, which is typically unobservable in macro- 
scopic samples. Advances in the synthesis and manipu- 
lation of biological matter have made it possible to study 
the response of individual molecules [if to mechanical 
pulling forces, also see Ref. [2J for a recent review. In 
these experiments, the external force is modulated follow- 
ing a predetermined protocol and the resulting work W 
exerted upon the system is measured. Consequently, the 
fluctuating nature of the thermodynamic quantity 'work' 
is revealed, as repeated measurements yield different re- 
sults depending on the microscopic trajectory traversed 
by the molecule. Similar studies have been performed on 
colloidal particles Q and mechanical oscillators H . Fluc- 
tuation theorems such as the Jarzynski equality [5] enable 
one to infer path-independent equilibrium free energy dif- 
ferences from the statistics of irreversible measurement 
protocols. The Jarzynski equality (JE) states that 



-W/k-eT 



-AF/k B T 



(1) 



where the ensemble average is taken over repetitions of 
the pulling protocol starting from an equilibrium state at 
temperature T, and AF is the difference in free energy 
for final and initial values of the external control param- 
eters. As we demonstrate experimentally, this equation 
applies also to a single-electron device operated in a di- 
lution refrigerator at temperatures of about 200 mK and 
below. An electrostatic gate drive plays the role of the 
time-dependent mechanical force that has been missing 
in previous studies of fluctuation theorems in electronic 
systems @ , and is crucial for studying JE ((T|) and Crooks 
fluctuation relations discussed below. 

The system that is the subject of our study is a single- 
electron box [?J as depicted in the scanning electron mi- 
crograph of Fig. 1(a). The box is formed by two metallic 
electrodes that are electrically connected through a tun- 
nel junction with a small electric capacitance Cj. The 
single degree of freedom that is of interest in the present 
experiment is the number n of electrons transferred from 



the left to the right electrode. The relevant n-dependent 
part [8| of the electrostatic energy of the idealized capac- 
itance system of Fig. 1(b) is then 



U = E c (n 2 - 2nn g ), 



(2) 



where E c — e 2 /(2Cs) is the characteristic unit of charg- 
ing energy of the box, e is electron charge, Ce = Cj + C g 

is the total capacitance of the box, C g — (C^ 1 + C^ 1 ) 
is the effective gate capacitance, and n g — C g V g /e is the 
gate charge in units of e. The charge number n changes 
by ±1 in the process of electron tunneling across the junc- 
tion. The tunneling is essentially instantaneous on the 
other relevant timescales in the problem j^j. In the tem- 
perature range of strong Coulomb blockade T <C E c /k&, 
thermal excitations of n are exponentially unlikely when 
the gate charge is an integer multiple of e relative to the 
background charge level of the box [tJ . Thus charge num- 
ber n can be controlled by the macroscopic gate voltage 
V g as demonstrated in Fig. 1(c). 

A typical setup for experiments in single-molecule sta- 
tistical mechanics is based on a laser-optical trap, in 
which the work W done on the molecule can be de- 
termined as an integral W — J. f FdX, where F is the 
measured force and A is a suitably-defined displacement 
variable. In contrast, the present experiment employs a 
calorimetric readout that yields the heat Q — W — AU 
deposited into the reservoir electrodes, where AU is the 
change in the internal energy of the system. In gen- 
eral, the relationship between Q and the dissipated work 
Wdis = W — AF appearing in the fluctuation relations 
depends on the charge state n of the system at the be- 
ginning and end of the gate protocol [8|, as illustrated 
graphically in Fig. 1(d). However, in the case of strong 
Coulomb blockade, the system can be made to virtually 
always follow a trajectory from n = to n = 1, in which 
case Wdi s = Q [Icj . Thus the Jarzynski equality assumes 
the form originally introduced in 11 | 



-Q/fesT 



(3) 



Importantly, the above formulation of JE can be applied 
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FIG. 1. (a) Scanning electron micrograph of the active 
area of the measured sample which shows metallic films fab- 
ricated on an oxidized silicon wafer by e-beam lithography 
and shadow evaporation technique 18]. Two shifted copies 
of the original resist mask pattern lie on the surface: Copper 
layer appears brighter compared to oxidized aluminum. Tun- 
nel junctions are formed in the overlap regions between the 
two films. The single-electron box is located on the left, and 
the SET electrometer at the top. The tips of two gate elec- 
trodes that are used to control the electrostatics of the box 
and the electrometer are visible at the left and right edges, 
(b) Simplified circuit diagram of the system. The isolated 
single-electron box is connected capacitively to its environ- 
ment through Cl and Cr and to the electrometer, illustrated 
in gray, (c) Experimental electrometer trace (black) together 
with the sinusoidal drive signal (red) applied to the control 
gate. The drive frequency is 1 Hz, and the amplitude is equal 
to one gate modulation period of the box. Close-up of the 
gray segment shows stochastic switching events as the drive 
signal passes the degeneracy point, (d) Energy level diagram 
of the system for the two lowest-energy charge states. Black 
parabolas represent the charging energy of the system in the 
states n = and n — 1 as a function of the extrnally controlled 
gate charge n g . Possible starting and ending configurations 
are indicated by open and filled circles, respectively, and the 
red annotation text indicates the dissipated work W — AF 
for each trajcetory type, illustrating that work separates into 
dissipated heat Q and change of internal energy. 



without detailed knowledge of the internal dynamics of 
the system, thus retaining its universality and usefulness. 
In this setting, one can furthermore write the Crooks 
fluctuation theorem [l2[ in terms of generated heat as 



P f{-Q) = -Q/k B T 

Pr(Q) 



(4) 



when the system is driven in forward (F) or reverse (R) 
directions, respectively. In our experiment, the calorime- 
try is performed by analyzing the signal from a capaci- 
tively coupled single-electron transistor (SET) electrom- 
eter in a manner outlined below [loj ]. The heat Q is 
obtained in units of the charging energy E c of the sys- 
tem. One can thus utilize the experimental distributions 
P{Q) in two ways: Using values of E c and T determined 
by independent means, validity of Eqs. © and (g]) can 
be tested. On the other hand, accepting Eqs. ([3]) or (jU, 
one can determine the ratio E c /k-QT, and furthermore 
find E c by multiplying this ratio with the independently 
measured temperature of the sample holder. 

Charge tunneling is governed thermally by the heat 
bath of phonons that couple to the excitations of con- 
duction electrons in the box electrodes. Hence, switching 
dynamics between the different charge states n is dissipa- 
tive. A quasiparticle with energy E tunneling from the 
left electrode to the right will extract an amount of energy 
E — hl from the left and deposit E — fiR to the right elec- 
trode. Hence, the total amount of energy deposited into 
the two electrodes equals the difference of the chemical 
potentials of the two electrodes, — "fl- The chemical 
potential difference is given by the change in energy U of 
Eq. @ due to change An = ±1 of the charge number n. 
We study the heat Q dissipated in the electrodes of the 
box in a gate control protocol n g (t) that ramps the gate 
voltage at an adjustable rate between the values n g = 
and n g = 1. At temperature T <C E c /ks, the box be- 
haves as a two-state system as it is driven between the 
states n — and n = 1. The n-trajectory consists of a 
random number N of successive back-and-forth tunnel- 
ing events and one can write the heat generated in such 
a trajectory as [ll[ 
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(5) 



where Pp and Pr are the probability distributions of Q 



where is the stochastic time instant of the fcth tun- 
neling event, and the sign is given by the change of n in 
the event. Equation (|5|) can be used to determine Q in 
units of E c for an individual n g (t) sweep. Because of the 
intrinsic randomness of the tunneling events, the heat Q 
fluctuates from one gate voltage ramp to another. In the 
experiment, the system was driven with a sinusoidal ex- 
citation corresponding to n g (t) = ^ — \ cos (2nft) with 
frequencies / ranging from 1 to 20 Hz. These frequen- 
cies are sufficiently slow so that the charge state always 
assumes its minimum energy value at the turning points 
of the drive, when n = n g {t) =0 and 1, respectively. 
Hence, each half cycle from to 1, and similarly from 1 
to 0, can be considered an independent realization of the 
gate control protocol. 

The preceding discussion is independent of the details 
of the charge tunneling rates in the single-electron box. 
However, in our case it is possible to analyze the flue- 



tuation relations also from a microscopic point of view. 
Charge transport through a Cu/AlO^/Al normal metal- 
insulator-superconductor (NIS) tunnel junctions occurs 
via thermally activated (TA) le events described by the 
orthodox theory [§], provided that (i) the tunneling re- 
sistance i?T of the junction is high compared to the resis- 
tance quantum i?x — 25.8 k£l, (ii) quasiparticles in the 
electrodes obey an equilibrium thermal distribution, and 
(iii) coupling of stray microwaves to the junction has been 
prevented by appropriate shielding and filtering in the 
construction of the sample stage and signal lines. Real- 
ization of these conditions in NIS single-electron devices, 
including the back-action from the capacitively coupled 
elecctrometer, has been studied in detail in recent years 
[l3l 14 1 . Based on these studies, non-thermal charge 
transport is expected to be negligible at least above tem- 
peratures of 150 mK. In particular, overheating of the 
superconducting electrode in the present design is di- 
minished by the fact that quasiparticle excitations can 
relax to the overlapping normal metal through the ox- 
ide barrier. Observed stochastic switching of the sys- 
tem between charge states at a fixed value of gate charge 
near degeneracy can be directly fitted to the tunneling 
rates predicted by the orthodox theory [l(|. We extract 
values A = 218 ± 3 /xeV, E c /k B = 1.94 ± 0.05 K, and 
Rt = 100 ± 13 MSI for the superconducing gap parame- 
ter, charging energy of the box, and tunneling resistance 
of the box junction, respectively. For the detector SET, 
we obtain Rt = 0.63 MSI and A = 211 /ieV from a fit to 
the measured I-V characteristics. 

The experimentally obtained Q distributions for drive 
frequencies 1, 2 and 4 Hz are presented in Fig. 2(a). The 
distributions were measured at a bath temperature of 
214 mK where the thermally activated tunneling rate at 
degeneracy was 70 Hz, which is well within the detector 
bandwidth of about 1 kHz. In addition, we have simi- 
lar data but in smaller quantities for driving frequencies 
from 5 Hz to 20 Hz. At frequencies higher than this, the 
observed distributions deviate significantly from the the- 
oretical prediction due to systematic errors arising from 
finite readout bandwidth and uncertainty in the event 
timing. On the other hand, driving frequencies lower 
than 1 Hz make the measurement susceptible to 1// type 
charge noise that is ubiquitous in metallic single-electron 
devices In Fig. 2(a), we show also the exact theoret- 
ical distributions based on charge kinetics described by 
the orthodox theory and sample parameters obtained in 
the manner described above. The theoretical and exper- 
imental distributions are in excellent agreement. 

To assess quantitatively the systematic error due to fi- 
nite detector bandwidth, we show in Fig. 2(a) also the 
distributions obtained from Monte Carlo simulations that 
incorporate a finite detector rise-time before threshold 
detection. Visually, the change in the shape of the dis- 
tribution functions appears small. Quantitatively, we 
can assess the accuracy of the readout by evaluating 
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FIG. 2. (a) Measured distribution of the generated heat 
at drive frequencies 1 Hz (black squares), 2 Hz (red circles), 
and 4 Hz (blue diamonds). The solid lines are exact theo- 
retical predictions for the independently determined sample 
parameter values. The dashed lines show the results of Monte 
Carlo simulations, where the finite bandwidth of the detec- 
tor was included in the model. Inset: P(Q)/P(—Q) ratio 
for the experimental distributions. Solid line shows the re- 
sult from the Crooks fluctuation theorem, (b), (c) First and 
second moments, respectively, of the Q distribution at dif- 
ferent drive frequencies and bath temperatures. Markers are 
experimental data, and solid lines are exact theoretical pre- 
dictions as in part (a). For the lowest bath temperatures, 
the theoretical curves have been calculated using a slightly 
elevated electron temperature to account for non-ideal ther- 
malization as discussed in the text. The temperature used 
in the calculation is given in brackets if it differs from the 
sample stage temperature. In panel (c), dashed lines repre- 
sent the distribution width inferred from the FDT formula 
{(Q ~ (Q)) 2 ) = 2fcT (Q) using the theoretical value of the 
first moment (Q). 



the exponential average (erQ' k ' aT \ , which equals 1 for 
the ideal thermally activated kinetics. From the Monte 
Carlo simulations, we obtain 1.006 for the 1 Hz and 2 Hz 
cases, and 1.012 for the 4 Hz case. For the experimental 
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distributions, evaluation of the Jarzynski average yields 
1.033 ± 0.003 (for 1 Hz drive), 1.032 ± 0.003 (2 Hz) and 
1.044 ± 0.004 (4 Hz), when using the independently de- 
termined E c /k,BT ratio as described above. The stated 
uncertainty is the unbiased estimate for the standard de- 
viation of the mean based on the observations, not in- 
cluding the uncertainty of the value of E c /k&T. As the 
relative uncertainty of the independent E c estimate is 
3%, the Jarzynski equality is shown to hold within ex- 
perimental accuracy, accounting for the 1% bias from fi- 
nite detector bandwidth. Conversely, starting from the 
assumption that JE holds for the experimental distribu- 
tions, we obtain an estimate E c /k-Q = 1.91 ± 0.03 K. 

The possiblity to evaluate numerically the theoretical 
Q distribution to a high accuracy enables us to assess the 
magnitude of sampling error in the experiment. From the 
acquired experimental data, the distribution could be de- 
termined for a range of Q/E c values where P(Q) > 10~ 3 . 
Evaluating numerically the contribution to JE outside 
this interval, we see that the deviation due to sampling 
errors is of the order 10~ 3 . 

Concerning the Crooks fluctuation theorem in the form 
of Eq. ([IJ , the forward and backward distributions coin- 
cide in the present case of a gate drive that is antisym- 
metric with respect to the degeneracy point. Hence, we 
present the experimental P(Q)/P(—Q) ratio as a func- 
tion of heat Q in the inset of Fig. 2(a). On a semilog- 
arithmic plot, one expects a linear dependence with the 
slope equal to l/kg,T, independent of frequency. The ex- 
perimental data adheres to this quite well, but there is a 
tendency towards less steep slopes for a lower frequency. 
This feature can be reproduced in our simulations by as- 
suming a broadening of the Q distributions due to an 
additive Gaussian noise having a r.m.s. amplitude of 
0.0035 E c = 0.032 k^T independent of driving frequency. 
We attribute the broadening to the residual background 
charge noise that induces fluctuations in the exact posi- 
tion of the degeneracy point. 

In Figs. 2(b), (c) we present the mean generated heat 
/i = (Q) and distribution width a = \J ((Q — /i) 2 ) for dif- 
ferent bath temperatures in the range 66 — 214 mK and 
drive frequencies 1 — 20 Hz. The experimental results 
agree well with the values obtained from exact theoret- 
ical distributions, which we calculated using the previ- 
ously determined sample parameters. The moments of 
Q distribution have been studied theoretically for a nor- 
mal state box in Ref. 11], and similar results hold for 



the first two moments in the present NIS case as well. In 
general, one expects Q = in the adiabatic limit / — > 0. 
The behavior is evident in the experimental data for the 
highest temperatures, whereas for low temperatures, the 
adiabatic threshold frequency is well below 1 Hz and thus 
inaccessible in the present experiment. Furthermore, for 
sufficiently low /, the box remains close to local equi- 
librium during the gate voltage drive. In this case, the 
transfer of heat Q into the box reservoirs can be viewed 



as a linear response to small deviations from the equi- 
librium caused by the drive. As a linear response, this 
process satisfies the classical fluctuation-dissipation the- 
orem (FDT) which takes the form a = y/2k-BTfi. This 
implies that in the adiabatic limit, the generated heat 
vanishes not only on average, but for all individual tun- 
neling trajectories. Local equilibrium also implies that 
the distribution of heat Q is Gaussian LLj similarly to 
all equilibrium thermodynamic fluctuations. In our ex- 
perimental results, this approximation holds only for the 
few lowest drive frequencies at 214 mK. Outside this re- 
gion, the distributions are strongly non-Gaussian. For 
the theoretical calculations at the lowest temperatures, 
we had to use somewhat higher temperatures than those 
indicated by the sample stage thermometer in order to 
reproduce the experimental data [see Fig. 2(b)]. At low 
temperatures, charge kinetics is expected to depart from 
the basic thermal activation model with a single heat 
bath, as discussed earlier. 

Thermodynamics of the driven transitions in electronic 
systems studied in this work will play an important role 
in the development of reversible information processing 
devices • We also envision this work to introduce a 
fruitful testbed of non-equilibrium fluctuation theorems 
in qualitatively new settings such as engineered environ- 
ments. It is possible to realize experimentally a regime 
where the charge transitions are dominated by coupling 
to an external nonequilibrium environment 13j . A fully 
superconducting box should enable the study of ther- 
modynamic fluctuations in true quantum regime, which 
is mostly an unexplored territory in experiments at the 
moment. 
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We monitor the time-resolved tunneling events by a capacitively coupled SET electrometer as in 



Refs. 



BE 



■3|. The electrometer is voltage biased and tuned to a charge-sensitive operation point. The 
current is read out with a room-temperature transimpedance amplifier. Cross-coupling between the 
gate of the electron box and the SET island is compensated so that the detector current assumes 
the same value for a fixed charge state n of the box regardless of the control gate voltage. At 
temperatures T < E c /k-Q, the detector current will display a characteristic Coulomb staircase as 
a function of V g . The conversion factor C g /e from the applied gate voltage V g to gate charge n g 
can be accurately determined from the observed period of the staircase. Slow drifts in background 
charge during data acquisition were compensated using a first order closed-loop control with a time 
constant of 100 s. 

To apply Eq. (5) of the main article, we have analyzed the experimental traces with a non- 
hysteretic threshold detector for determining the switching times . Use of the threshold detection 
rejects most of the charge noise in the readout, and the residual noise appears as small jitter in the 
switching times. Besides intrinsic charge noise, the finite bandwidth of the charge readout leads 
to two additional error sources. First, there is a constant time delay between the instant of the 
tunneling event and the time when the filtered detector signal crosses the threshold level. In our 
setup, most of the time delay and bandwidth limitation originates from the digital low-pass filter 
employed to reduce the current noise of the readout prior to the threshold detection. The constant 
time delay has been compensated for in the analysis of the experimental traces. Secondly, finite 
bandwidth can cause back-and-forth transitions occurring in a rapid succession to be missed alto- 
gether [J]. To assess the effect of the finite detector bandwidth on the signal processing algorithms, 
we performed Monte Carlo simulations using a detector step response function extracted from the 
experimental traces. The results of these simulations are discussed in the main article. 

We have measured the rate of thermally activated (TA) transitions as a function of gate charge 
and bath temperature. Part of that data is presented in Fig. SI (a). The TA one-electron tunneling 
rates obey the detailed balance, i.e., T(—E) = e~ E ^ kBT T(E), where E is the chemical potential 
difference between the electrodes. Hence, E c can be determined from the tunneling rates observed 
at a fixed gate charge n g by utilizing A;BTln(ro^i/ri^o) = 2E c (n g — 1/2). We have estimated 
E c = 1.94 ± 0.05 K in this manner from data obtained at different temperatures, see Fig. Sl(b). 
Furthermore, for A and k^T <C A, which is the case in the experiment, one can write the 

approximation 




FIG. SI. Observed rate of thermally activated tunneling as a function of the gate charge offset from 
degeneracy at bath temperatures 197, 218, and 238 mK. Open and filled symbols correspond to —> 1 and 
1 — >• transitions, respectively. Solid lines show the theoretical prediction according to Eq. (Q]) with fitted 
sample parameters, (b) Estimates for the box charging energy E c obtained by application of the detailed 
balance to the observed To-^i/ri^o ratio as a function of gate offset at different temperatures, (c) Observed 
tunneling rate at degeneracy as a function of temperature (square markers), and best fit to Eq. (flj (solid 



where A is the superconducting gap parameter, and Rt is the tunneling resistance. We have 
determined Rt and A by a nonlinear fit to the tunneling rates observed at different temperatures 
at degeneracy, i.e., E = above. The observed rates are plotted in Fig. Sl(c) together with the 
rates predicted by Eq. (1) for the fitted sample parameters. 

Finally, for a quantitative assessment of the Wdy = Q approximation employed in the interpre- 
tation of the data, we have calculated numerically the full Q distributions also for the less likely 
0—7-1,1—7-1, and 1—7-0 trajectory types using the experimental sample and operation parameters 
stated in the main article. Given that E c /ks = 1.94 K, the probability to start (finish) the gate 
trajectory in a state other than n = (n = 1) is approximately 10~ 4 at the highest temperature 
of 0.214 K for which we present experimental data. Hence, the distribtions for Wdis an d Q are 
expected to differ only at the 10 -4 level. The Q distributions can be combined into the full Wdis 
distribution using the relations presented in Ref. [5j. The resulting distributions are shown in 
Fig. S2, and one observes that the distributions for Wdis and Qo->i coincide, hence validating the 
approximation used in the article. 
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FIG. S2. Full simulated Q distributions for the four different trajectory types weighted with their occurrence 
probability plotted as dashed lines. The simulations were made for a / = 5 Hz sinusoidal drive at T = 
214 mK, and the sample parameters and operation conditions correspond to those stated in the main article. 
In Fig. 2(a) of the main article, only the Qo->i distributions are presented. Here, the full distribution of 
W—AF is overlaid on the graph with a solid line [with (W— AF)/E C as the x axis], and it is indistinguishable 
from that of Qo-vi- 
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